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Abstract 

Kramers-Grote-Hynes factor is the key nonequilibrium contribution to rate constant of a reaction 
over and above the transition state theory rate in the spatial limited regime. Wolynes in eighties 
introduced a quantum correction to the overall rate coefficient. This is responsible for tunneling and 
quantum enhancement of rate at low temperature. However, its validity is restricted to activated 
tunneling region or above crossover temperature. Based on a quantum formulation of the normal 
mode analysis, we suggest a generalization of Wolynes factor and a consequent multidimensional 
transition state rate expression which are valid in the deep tunneling region down to zero degree 
Kelvin. 
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I. INTRODUCTION 



Kramers- diffusion model of chemical reactions proposed in 1940 forms the dynamical 
basis of modern rate theory of activated processes. The seminal and essential content of 
this nonequilibrium formulation is the inclusion of dependence of rate constant of a reaction 
on viscosity or friction of the reaction medium. Based on the classical theory of Brownian 
motion in phase space Kramers derived the expressions for nonequilibrium steady state 
distribution functions to work out the rate coefficients in the two different limiting situations 
and showed that the rate varies linearly in the weak dissipation regime and inversely in 
the high dissipation regime with friction. With the advent of ultrafast lasers and time- 
resolved detection techniques since late seventies, experimental confirmation of Kramers' 
theory provided a new impetus for further development in chemical dynamics and condensed 
matter physics ^ 3 ' 4 ^ 6 . 

While Kramers' theory is based on classical Markovian description, eighties saw a num- 
ber of important generalization of Kramers theory in several directions by a number of 
group o 7 ! 8 ! 9 ! 10 ! 11 ! 12 ! 13 ! 14 ! 15 . An important endeavor in this context is the inclusion of a quan- 
tum equilibrium factor, known as Wolynes facto r 16 i 17 , in the generalized Kramers' or equiva- 
lently multidimensional transition state theory rate constant. This term allows tunneling as 
a possible decay route of the metastable state in addition to the usual thermally activated 
processes and is responsible for quantum enhancement of rate at low temperature. However, 
the validity of the Wolynes factor is restricted to the activated tunneling regime or above 
crossover temperature. Based on a positive definite Wigner function formalism within a nor- 
mal mode description we generalize Wolynes factor and thereby multidimensional transition 
state rate expression which are valid in the deep tunneling region and vacuum limit. 

The outlay of the paper is as follows. In the following Sec. II we discuss a quantum 
Langevin equation which allows us to realize a c-number Hamiltonian amenable to a normal 
mode analysis. The equilibrium theory in terms of Wigner distribution function to formulate 
a quantum counterpart of multidimensional transition state theory (TST) has been presented 
in Sec. III. We work out the rate expression and the Wolynes factor in the appropriate limits. 
The paper is concluded in Sec. IV. 
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II. QUANTUM LANGEVIN EQUATION 



We consider a particle of unit mass coupled to a medium comprised of a set of harmonic 
oscillators with frequency cjj. This is described by the following Hamiltonian: 



d 4 + m+ f{f + i(. 



UiXi q 

U)i. 



(2.1) 



Here q and p are co-ordinate and momentum operators of the particle and the set {£i,pi} 
is the set of co-ordinate and momentum operators for the reservoir oscillators coupled linearly 
to the system through their coupling coefficients Cj. The potential V(q) is due to the external 
force field for the Brownian particle. The co-ordinate and momentum operators follow the 
usual commutation relations [q,p] — ih and [£i,Pj] = ihSij. Here masses have been assumed 
to be unity. 

Eliminating the reservoir degrees of freedom in the usual way we obtain the operator 
Langevin equation for the particle, 



lit) + / df 7 (t - i(f) + V(q) = F(t) 



(2.2) 



where the noise operator F(t) and the memory kernel j(t) are given by 



—Xj(0) — q(0) > cosojjt + —pj(0) sinujjt 



(2.3) 



and 



N a 
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(2.4) 



The Eq. fj2.2)l is the well known exact quantized operator Langevin equation for which the 
noise properties of F(t) can be derived by using a suitable initial canonical distribution of 
the bath co-ordinate and momentum operators at t = as follows; 



(F(t)) QS = (2.5) 
^{(Ht)F(t'))QS+(F(t')F(t)) QS } = coth ^B 008 Wi (t-f) (2.6) 

9=1 ■? 
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where (..-)qs refers to quantum statistical average on bath degrees of freedom and is 
defined as 



(d) QS = ™^P(zgWMj (2.7) 
Trexp(- H hath /k B T) 

for any operator 0({(u?/cj)xj - q}, {pj}) where H h&th (Y!i=i\Pi V 2 + V 2 ^^ _ f~"<?) 2 ]) at 
t = 0. By Trace we mean the usual quantum statistical average. Eq. (j2.6J) is the fluctuation- 
dissipation relation with the noise operators ordered appropriately in the quantum mechan- 
ical sense. 

To construct a c-number Langevin equatio n 20 ' 21 we proceed from Eq. (j2.2j) . We carry out 
a quantum mechanical average of Eq. ()2.2|) 

(q(t)} + fdt'^t-t'Mt')} + (V'(q)) = (F(t)) (2.8) 
Jo 

where the quantum mechanical average (...) is taken over the initial product separable 
quantum states of the particle and the bath oscillators at t — 0, |0){|ai)|a!2) . . . |ajv)}- 
Here \<p) denotes any arbitrary initial state of the particle and \cii) corresponds to 
the initial coherent state of the z-th bath oscillator. \oti) is given by |a;j) = 
exp(— |aj| 2 /2) Y^ i= o( a T I V^-)\ n i) i a i being expressed in terms of the mean values of 
the shifted co-ordinate and momentum of the i-th oscillator, {(a; 2 /c i )(x i (0)) — (q(0))} = 
(y/h/2ui)(ai + a*) and (p«(0)) = i^/huJi/2(a* — «j), respectively. It is important to note 
that (F(t)) of Eq. (j2.8)l is a classical-like noise term which, in general, is a non-zero number 
because of the quantum mechanical averaging and is given by ((F(t)) = f(t)); 



m = E — fo(°)> - w)> 4 cos ^ + -<Pi(o)> dwjt 

It is convenient to rewrite the c-number equation (|2.8|) as follows; 



(2.9) 



($(*)> + f dt'^t-t'mt')) + (V'(q)) = f{t) (2.10) 
Jo 

To realize f(t) as an effective c-number noise we now introduce the ansatz that the 
momenta (pj(0)) and the shifted co-ordinates {(ix> 2 /cj)(xj(0)) — (q(0))}, {pj} of the bath 
oscillators are distributed according to a canonical distribution of Gaussian form as 
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fe(o)) 2 + |{^(%(o))-(g(Q))} 2 



(2-11) 



so that for any function of the quantum mechanical mean values Oj{(pj(0)) , 
((u)j/cj)(xj(0)) — (q(0)))} of the bath the statistical average (. . .)s is 



(Oj)s 



O j Pj c%(0)> d{{uy Cj ){xM) - (?(0)» 



(2.12) 



Here fij indicates the average thermal photon number of the j-th oscillator at temperature 
T and fij(ujj) = l/[exp {fuVj/ksT) — 1] and N is the normalization constant. 

The distribution ()2.11|) and the definition of statistical average ()2.12|) imply that f(t) 
must satisfy 



That is, c-number noise f(t) is such that it is zero-centered and satisfies the standard 
fluctuation-dissipation relation (FDR) as expressed in Eq. fl2.6|) . It is important to emphasize 
that the ansatz (|2.11|) is a canonical Wigner distribution for a shifted harmonic oscillato r 18 ! 19 
and an exact solution of Wigner equation, which remains always a positive definite function. 
A special advantage of using this distribution is that it remains valid as pure state non- 
singular distribution function at T = 0. Furthermore, this procedure allows us to bypass 
the operator ordering prescription of Eq. ()2.6|) for deriving the noise properties of the bath 
in terms of fluctuation-dissipation relation and to identify f(t) as a classical looking noise 
with quantum mechanical content. 

We now return to Eq. (j2.1U|) to add the force term V'((q)) on both sides of Eq. (j2.1U|) and 
rearrange it to obtai n 20 ! 21 



(f(t))s = 



(2.13) 



and 




(2.14) 



q = p 





(2.16) 
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where we put (q(t)) = q(t) and (p(t) = p{t) for notational convenience and 

Q = V'((q)) - (V'(q)) (2.17) 

represents the quantum correction due to the system degrees of freedom. Eq. (|2.16[) offers 
a simple interpretation. This implies that the quantum Langevin equation is governed 
by a c- number quantum noise f(t) originating from the heat bath characterized by the 
properties ()2.13|) and ()2.14|) and a quantum fluctuation term Q 20 i 21 i 22 i 23 characteristic of the 
non-linearity of the potential 2 ^. 

Referring to the quantum nature of the system in the Heisenberg picture, we write 22,23 . 



q(t) = q + 5q (2.18) 
p(t) = p + 5p (2.19) 

where (q)(= q) and (p){= p) are the quantum-mechanical averages and 5q, 5p are the 
operators. By construction (5q) and (5p) are zero and [6q, 5p] = ih. Using Eqs. (j2.18|) and 
(|2.19|) in (V'(q)) and a Taylor series expansion around (q) it is possible to express Q as 

Q(Q, (6q n )) = - E ^ n+1 \q){Sr) (2.20) 

n>2 

Here V^ a \q) is the n-th derivative of the potential V(q). For example, the second order 
Q is given by Q = —\V"(q){5q 2 ). The calculation of Q 2 fli 21 i 22 i 23 therefore rests on quantum 
correction terms, (Sq n ) which are determined by solving the quantum correction equations 
as discussed in the following subsection. 

A. Quantum correction equations 

We now return to operator equation ()2.2|) and put ()2.18|) and ()2.19|) and use of Eq. (j2.10|) 
to obtain the following operator equations 

5q = 5p (2.21) 
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sp+ [\(t - 1) 5p(f) at + v"( q )5q+j2^y (n+1 \i) m n -(sq n )) 

J n >2 n - 

= F{t) - /(*) ( 2 -22) 

Eqs. (2.21-2.22) form the key element for calculation of quantum mechanical correction 
(J2.20)) due to nonlinearity of the system potential. Depending on nonlinearity of the potential 
and memory kernel we consider the following cases separately. 

1. Arbitrary memory kernel and harmonic potential 

We consider the friction kernel to be arbitrary but decaying and the potential as harmonic 
for which the derivatives of V(q) higher than second vanishes. After quantum mechanical 
averaging over the initial bath states J] t ~i {|«i(0))} only, the Eas. (l2.21H2.2^|) yield 

5q = 5p (2.23) 

5p = - I j(t- t') 5p(t') dt' - u 2 5q (2.24) 
Jo 

where u 2 = V"(q), a constant and u is the frequency of the harmonic oscillator. The 
above operator equation can be solved exactly by Laplace transform technique to obtain 

Sq(t) = 5q{0) (l-uj 2 [ CJt')dA + 5p(0) CJt) (2.25) 



where C p (t) is obtained as the inverse Laplace transform of C p (s) given by 



°p( s ) = 2^ ^TTT 2 (2-26) 

with j(s) as 



7(s) = / -f(t)e- st dt (2.27) 
Jo 

is the Laplace transform of the friction kernel 7(i). After squaring and quantum mechan- 
ical averaging over arbitrary initial system state \<p), Eq. (j2.25|) yields 
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(8q 2 (t)) = (6q\0) > C 2 (t) + (8p 2 (0) > C 2 p (t) 

+ C p (t) C q (t) {(5p(0) 5q(0)) + (Sq(0) mm (2.28) 

with C q (t) = l-u 2 f* C p (t')dt'. 

2. Exponential memory kernel and arbitrary potential 

This corresponds to a commonly occurring situation for which the bath modes are as- 
sumed to follow a Lorentzian distribution characterized by a density function p(u>) such that 
for Eq. (j2.4j) in the continuum limit we write 

c 2 M , s 2 r , 

where V is the dissipation constant in the Markovian limit and r c refers to correlation 
time. The memory kernel is exponential and 7(t) is given by 

7 (t) = - e -l*l/^ (2.30) 

In this case the quantum correction equations (|2.21j) and ()2.22|) after quantum mechanical 
averaging over bath states may be rewritten 

5q = 5p (2.31) 



5p = -V"(q) 6q - — V {n+l \q) (6q n - (5q n )) + 5z (2.32) 

n>2 U ' 

r i 

5z = 5p 5z (2.33) 

where we have introduced an auxiliary operator 5z to bypass the convolution integral in 
Eq. (j2.22|) . Making use of the Eqs. (j2.31H2.33|) we derive the quantum correction equations 
upto second order to obtain 



. 2 

{8q ) = (5q5p + 5p5q) 
8 



(2.34) 



(6qSp + 5p5q ) = 2 (5p 2 ) - 2 V"{q)(5q 2 ) + (5q5z + 5z5q) - V"'{q)(Sq 3 ) 



(2.35) 



(5p 2 ) = -V"(q)(5q5p + 5p5q) + (5p5z + 5z5p) - -V"'(q)(5p6q 2 + 5q 2 5p) (2.36) 

(6p6z + 6z5p ) = - V"(q)(5q5z + 5z5q) - ^V"\q){5q 2 5z + 5z5q 2 ) 

9F 1 

(Sp 2 ) (5p5z + Sz5p) + 2 (Sz 2 ) (2.37) 

Tc T c 

r i 

(5q5z + 5z5q ) = (5p5z + 5z5p) (5q5p + 5p5q) (5q5z + 5z5q) (2.38) 

T c T c 

(Sz 2 ) = — (5pSz + 5z5p) - — (5z 2 ) (2.39) 

T c T c 

Discarding the third and higher order terms in Eqs. (j2.34M2.39|) we obtain a set of closed 
equations which can be solved numerically along with Eqs. ()2.15j) and (|2.16jl . using suitable 
initial conditions for calculating the leading order (second order) contribution to quantum 
corrections in Q.A standard choice of initial condition a 22 i 23 corresponding to minimum un- 
certainty states is (5q 2 (0)) = h/2u, (5p 2 (0)} = hw/2 and (5q(0)6p(0) + 5p(0)5q(0)) = h with 
the other moments fEqs. (j2.37M2.39|0 being set at zero. The procedure may be extended 
to include higher order quantum effects without difficulty and may be easily adopted for 
numerical simulation of quantum Brownian motion^. 

B. Calculation of quantum statistical averages 

Summarizing the discussions of the last two sections A and B we now see that the 
Langevin dynamics in c-numbers can be calculated for a stochastic process by solving 
Eqs. fl2.15|) and 1)2.16)1 for quantum mechanical mean values simultaneously with quantum 
correction equations which describe quantum fluctuation around these mean values. In prin- 
ciple for nonlinear system the equations for quantum corrections constitute an infinite set of 
hierarchy which must be truncated after a desired order, in practice, to make the system of 
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equations closed. Care must be taken to distinguish among the three averages, the quantum 
mechanical mean (0)(— O), statistical average over quantum mechanical mean (O)s and 
the usual quantum statistical average (O)qs as discussed in Sec. II. To illustrate the relation 
among them let us calculate, for example, the quantum statistical averages (<z)qs, (q 2 )Qs 
and (q 2 p) QS - By (l2~THj) and (f2~TTJl) we write 



q = q + Sq 
>QS = (q + Sq) QS 

= (q)s + ({Sq)) s = (q)s 



(2.40) 
(2.41) 



Again 



(q 2 ) Q s = ((q + Sq) 2 ) QS 

= (q 2 )s+((5q 2 ))s (2.42) 

In the case of harmonic potential, (Sq 2 ) as given by (|2.28|) is independent of q or p so 
that one may simplify ()2.42|) further as 

(q 2 )QS=(q 2 )s+(5q 2 ) (2.43) 

In Ref.— the explicit exact expressions for (c[ 2 )qs and (p 2 )qs have been derived for har- 
monic oscillator and they are found to be in exact agreement with those of Grabert et alM. 
For anharmonic potential, however, one must have to use ()2.42|) to carry out further the 
statistical average over (Sq 2 ) i.e ((Sq 2 ))s, since it is a function of stochastic variables q and 
p according to quantum correction equations. Furthermore we consider (q 2 p)Qs 



(q 2 p)qs = ((<? + Sqf(p + 5p)) QS 

= (q 2 p)s + (P (Sq 2 ))s + ((5q 2 5p)) s + 2(q(6q 6p)) s (2.44) 

The essential element of the present approach is thus expressing the quantum statistical 
average as the sum of statistical averages of set of functions of quantum mechanical mean 
values and dispersions. Langevin dynamics being coupled to quantum correction equations, 
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the quantum mechanical mean values as well as the dispersions are computed simultaneously 
for each realization of the stochastic path«&. A statistical average implies the averaging over 
many such paths similar to what is done to calculate statistical averaging by solving classical 
Langevin equation. Before leaving this section we mention a few pertinent points. 

First, the distinction between the ensemble averaging by the present procedure and by 
the standard approach using Wigner function is now clear. From Eq. (j2.44j) we note that, 
for example, 



where W(q,p) is the Wigner function for system operators. (This is not to be confused 
with the Wigner function we introduced in (j2.11)) for the bath oscillators.) 

Second, our formulation of the Langevin equation coupled to quantum correction equa- 
tions belongs to quantum stochastic process derived by c-number noise, which is classical-like 
in form. Its numerical solutions can be obtained in the same way as one proceeds in a clas- 
sical theory 20 . 

Third, quantum nature of the dynamics appears in two different ways. The heat bath 
is quantum mechanical in character whose noise properties are expressed through quantum 
fluctuation-dissipation relation. The nonlinearity of the system potential, on the other hand, 
gives rise to quantum correction terms. Thus the classical Langevin equation can be easily 
recovered (i) in the limit hio <C fcsT to be applied in the Eq. 1)2.14)1 so that one obtains the 
classical fluctuation-dissipation relation and (ii) if the quantum dispersion term Q vanishes. 

C. c-number Hamiltonian and normal mode description 

The c-number Hamiltonian corresponding to Langevin equation ([2.15)1 and ([2.16)1 is given 



Note that the above Hamiltonian is different from our starting Hamiltonian operator 
(12.1)1 because of the c-number nature of ()2.46j) . {xi,pi} are the quantum mean value of the 
co-ordinate and the momentum operators of the bath oscillators. 




(2.45) 



by 
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The spectral density function is defined as 



N 2 



(2.47) 



We now assume that at g = 0, the potential V(q) has a barrier with height such that 
a harmonic approximation around q = leads to 



7(g) = V x - -u 2 b q 2 + V 2 (q) 



(2.48) 



where 



V"(q) | g= o, refers to the second derivative of the potential V(q). u b is the 
frequency at the barrier top and V 2 (q) is the non-linear part of the classical potential and 
is given by V 2 = ^ n>3 ^ 9 Q^f 1 | 9 =o Q n - With Eq. (|2.48j ) the quantum correction part in the 
Hamiltonian Eq. (|2.4fij) becomes 



^V^(g)(5g") = -^B 2 + V 3 (q) 



(2.49) 



n>2 



where B n = (5q n ); V3 = ^2 n>2 9 ■ Note that we have introduced an approximation 
by putting a bar over quantum dispersion (5q n ) to indicate its time average since it is 
sufficient to consider the energy loss of the system mode averaged over one round trip time, 
i.e., the time required to traverse from one turning point of the potential well to another 
and back. Putting (|2.48j) and (|2.49|) in the Hamiltonian (|2.46j) we obtain 



H = H + V N (q) 



(2.50) 



where we have decomposed the Hamiltonian in the harmonic part H and the anharmonic 
part Vjv(g) as 



+ 



UiXi q 



(2.51) 



and 



V N (q) 
and V± 



V 2 (q) + V 3 (q) 

B 2 n 
T 



(2.52) 



V x - —u 
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and V%(q) are therefore classical and quantum anharmonic contributions to total 
anharmonic part of the Hamiltonian. The separability of the c-number Hamiltonian in the 
quadratic and nonlinear parts allows us to make a normal mode transformation to convert 
the quadratic Hamiltonian into a diagonal form. The method of normal mode analysis has 
been used extensively by Pollak and co-workers^. 

Following Pollak, we diagonalize the force constant matrix T of the Hamiltonian ()2.51|) 
with the matrix U 

UT = X 2 U (2.53) 
where U provides the transformation from old co-ordinates to the normal co-ordinates 



/ a \ 



P 
Vi 

J/2 



/ a \ 



U 



Q 

x 2 



(2.54) 



\y N J 



\ x N ) 



The c-number Hamiltonian of the unstable normal co-ordinate is given by 



Ho = \f> 2 + Vf - \\y + £ \ [vl + XM) (2-55) 

i=l 

The eigenvalues and are expressible in terms of the coupling constant of the system 
and the bath implicitly as follows: 



N 



c 2 



X 



N 



1,2...N 



(2.56) 
(2.57) 



where (J2.56j) and (|2.57|) correspond to normal mode frequencies of the unstable mode and 
the i-th bath oscillator respectively. 
The transformation ()2.54|) implies 
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N 

q = MOO p + ^ U j0 yj 
3=1 



[2.58) 



and it has been shown 26 that woo and Ujo may be expressed as 



u 



00 



A' 



3=1 



'rf + A 2 ) 2 



(2.59) 



and 



u 



jo 



N 



-, -1 



(2.60) 



Making use of the spectral density function 1)2.47)1 and Laplace transformation of ~f(t), 
Eq. (J2.56j) and Eq. 1)2.59)1 may be written in the continuum limit as 



l + 7(A fe )/A 6 



(2.61) 



and 



u, 



00 



du- 



uj J{uj) 



-I -1 



(2.62) 



Wo (Ag + ^ 2 ) 2 . 

The two important identities in relation to orthogonal transformation matrices and the 
associated frequencies may be noted here for the dynamics at the barrier top and at the 
bottom of the well; 



N 



N 



U), 



IK = A >II A ? 



(2.63) 



i=l 



i=l 



and 



N 



N 



n^ 2 = A ori A ' 



(2.64) 



i=l 



i=l 



Here ujq and Ao are the frequencies of the system at the bottom of the well in the original 
co-ordinate and normal co-ordinate respectively. Similarly A* corresponds to the normal 
mode frequencies of the bath oscillators corresponding to a normal mode Hamiltonian at 
the bottom of the well, 
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ttI 1 -/2 i 1 \2 /2 i \ 1 -12 i 1*2 /2 

H o = 2 p + 2 oP + I 2 % 2 4 % 

U=i 



^2.65) 



Here // and p' are coordinate and momentum of system mode respectively and y[ and y\ 
are coordinate and momentum of i th bath oscillator respectively at the bottom of the well 
in the normal coordinates. Aq and A; are given by 2 § 



X 2 - , , 2 
A — LU 



N 



-, -1 



N 



C 2 



, i = l,2..JV 



(2.66) 
(2.67) 



III. C-NUMBER QUANTUM VERSION OF MULTIDIMENSIONAL TST 

To begin with we consider the particle to be trapped in a well described by a potential 
V(q). In the normal mode description of (N + 1) oscillators according to the Hamilto- 
nian (J2.55|) the bath modes and the system mode are uncoupled. Considering the un- 
stable reaction co-ordinate to be thermalized according to the Wigner thermal canonical 
distributio n 1 ^ 1 ? of N uncoupled harmonic oscillators plus one inverted we have 



P 



eq 



cxp 



i I/t 1 \2_2 

hX in (X ) + |) 



N 



H exp 



i=l 



\v1 + \^yf 



2 it>i 



hX l (n l (X i ) + |) 



(3.1) 



z is the normalization constant. As usual this can be calculated using the distribution 
function inside the reactant well. For this it is necessary to consider the normal mode 
Hamiltonian at the bottom of the well expressed as H' Q in the Eq. ()2.65|) . The corresponding 
distribution in the well is 



eq 



exp 



oP + 2 A 0P 



HX (n (X ) + \) 
which can be normalized to obtain 



n 

i=i 



exp 



(3.2) 



N 



A; 



27r^Ao(rr (Ao) + \) f = \ 2nM. i (n i (A i ) + §) 
The identity relation (|2.64j) can be used to transform to the following form 



(3.3) 
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U 



N 



to. 



2irh\ (n (\ ) + \) f \ 2nhh i {n i {k i ) + f) 
Putting Eq. ()3.4|) in Eq. ()3.1|) we obtain after integration over the stable modes 



P, 



eq 



2irh\ (n (\ c 



where 



TT — X exp 



1 n 2 4- \A 1 \ 2 n 2 

ftA (n (A ) + |) 



(3.4) 



(3.5) 



The total energy of the unstable mode is 



(3.6) 



E 



(3.7) 



2' 1 2 

The prime quantity for determination of rate constant is the distribution of energy of 
the unstable mode. Thus going over to an energy space so that the co-ordinate p, p are 
transformed to t, E, respectively with unit Jacobian, the equilibrium distribution function 
(|3.5|) is given by, 



fe q (E) 



A,, 



TT — X exp 



E 



(3.8) 



2irh\ (n (\ ) + 3) ^> 

The above distribution is valid for the energy of the unstable mode E > as well as 
E < V*. 



A. Quantum multidimensional TST rate 

The rate of activated barrier crossing in terms of the equilibrium probability becomes 

f{E) dE (3.9) 

vt 

As the unstable mode remains uncoupled from the stable modes the former mode behaves 
deterministically and the recrossing does not occur in this case. 

Making use of the distribution ()3.8|) in (J3.9|) we obtain the rate constant 
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The above expression corresponds to a quantum multidimensional transition state rate 
constant. This is central result of this section. Apart from usual Kramers-Grote-Hynes 
term \b/ujb and u /27i, the term arising out of classical transition state result, it contains 
two important factors. First, an exponential Arrhenius term where the usual thermal factor 
ksT is replaced by ^A (n (Ao) + \) includes quantum effects due to heat bath at very low 
temperature. In the high temperature limit it reduces to ksT and one recovers the usual 
Boltzmann factor. This term is essentially an offshoot of a description of thermal equilibrium 
by a canonical Wigner distribution of harmonic oscillator heat bath. Second term x can 
be identified as the quantum correction to Grote-Hynes factor or more precisely a vacuum 
corrected generalized Wolynes contribution for quantum transmission and reflection for the 
finite barrier. While usual Wolynes term takes care of the quantum effects at the higher 
temperature the factor \ incorporates quantum effects at arbitrary low temperature. In what 
follows we show that the usual Wolynes term and the quantum rate valid above cross-over 
temperature can easily be recovered from x in the appropriate limit. 

B. Derivation of Wolynes factor from \ and limiting rate expression at low tem- 
perature above cross-over 

We begin by noting that n(x) in x which is given by 



is the Bose distribution n (x) = ( e hx / kBT -I)' 1 . Neglecting the vacuum contribution 1/2 




(3.11) 



from the terms like hx(n(x) + |) and keeping only the leading order quantum contribution 
we obtain 




(3.12) 



Therefore x reduces to S (say) the Wolynes factor 
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From identities (gUSD and (glElD it follows 



N 



iV 



1=1 

JV 



Aft 



i=i 

iV 



1=1 



A„ 



(3.14) 
(3.15) 



i=i 



respectively, and we have 



N 

n 



Aj _ Ao^ 



(3.16) 



Making use of the relation ()3.16|) in ()3.13|) we obtain 



N 



Ao^6 tt sinh(^,Aj/2A; B T) 
c^oAfe smh(h\i/2kBT) 



(3.17) 



Furthermore (A ^ b )/(u;oAfc) can be rewritten as ^ r^T^f^] which may be approximated 
in the form [uji/ujq) sinh(^A°/2fcaT) • Eq. (|3.17j ) then reduces to 



A' 



„ Ub sinh(/lAo/2A;£T) -j-j- sinh(^Aj/2A;£T) 



(3.18) 



ujq smh(h\b/2kBT) -M- sinh(/iAj/2/csT) 

This is the wellknown Wolyne a 1 ^ 1 ? expression derived in eighties as a higher temperature 
equilibrium quantum correction to Kramers-Grote-Hynes dynamical factor to Kramers' rate. 
Both x an d the Wolynes factor become unity in the classical limit. We conclude by noting 
that unlike Wolynes factor H, \ is valid below cross-over temperature. The quantity 5 ap- 
proaches unity when the temperature far exceeds the cross-over temperature. The Wolynes 
factor may be approximated to yield a form at low temperature (valid above cross-over 



temperature) as exp 



h 2 K+^q) 

24 (fc s T) 2 



This term is responsible for well known T 2 enhancement 



of rate due to quantum effects. With this form of 5 and exp 



reducing to 



exp 
earlie^ 



k B T 



h\ {n Q (\ )+\) _ 

above the cross-over regime we obtain from ()3.10|) the rate expression derived 
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2n u b 



exp 



ft H + ^o 2 ) 

24 {k B TY 





r i 


exp 


i k B T\ 



(3.19) 



Our general expression ()3.10j) valid in the deep tunneling regime, i.e., down to vacuum 



regime. 



IV. SUMMARY AND CONCLUSION 



Based on a quantum Langevin equation we have constructed a c-number Hamiltonian for 
a system plus N-oscillator bath model. This allows us to formulate a normal mode analysis 
to realize a c-number version of the multidimensional transition state theory and to derive 
a quantum expression for the total decay rate of metastable state. The result is valid for 
arbitrary noise correlation and temperature down to vacuum limit. The following pertinent 
points are noteworthy. 

(i) We have shown that the expression for quantum rate coefficient is a product of four 
terms , e.g. , classical well frequency, Kramers'-Grote-Hynes factor, a vacuum corrected 
or generalized Wolynes factor representing quantum transmission and reflection, an expo- 
nential term corresponding to Wigner canonical thermal distribution i.e. the generalized 
Arrhenius term. Of these the Wigner term and the vacuum corrected Wolynes term refer 
to equilibriation of quantum particle in the well and therefore corresponds to a c-number 
multidimensional transition state result. Since the distribution unlike the Boltzmann is valid 
even as T — ► 0, the quantum effect due to heat bath can be well accounted by these terms 
even below the activated tunneling regime. 

(ii) The classical limit of the quantum rate expression depends on Wigner, generalized 
Wolynes factor. It is easy to see that they reduce to Arrhenius factor, unity respectively in 
the limit tko <C ksT. 

(iii) The present theory takes care of activation and tunneling within a unified description 
and is equipped to deal with the rate at temperature down to vacuum limit. This is a distinct 
advantage over path integral Monte Carlo method since numerically the relevant propagator 
poses serious problem as the temperature approaches absolute zero. 

The quantum theory as presented here is based on a canonical quantization procedure 
and a description of equilibrium by positive definite Wigner's thermal distribution for the 
bath rather than path integral or master equation formalism. The approach has been used 
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and are in use in many related issues^. The systematic improvement can be made by taking 
care of the quantum corrections of higher orders. 
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